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ABSTRACT 

While BL Lacertae objects are widely believed to be highly beamed, low-luminosity radio galaxies, 
many radio-selected BL Lacs have extended radio power levels and optical emission lines that are too 
luminous to be low-luminosity radio galaxies. Also, Stocke & Rector discovered an excess of Mgll 
absorption systems along BL Lac siglitlines compared to quasars, suggesting that gravitational lensing 
may be another means of creating the BL Lac phenomenon in some cases. We present a search for 
gravitationally-lensed BL Lacs with deep, high-resolution, two-frequency VLA radio maps of seven 
lensing candidates from the 1 Jansky BL Lac sample. We find that none of these objects are resolved 
into an Einstein ring like B 0218-1-357, nor do any show multiple images of the core. All of the lensing 
candidates that were resolved show a flat-spectrum core and very unusual, steep-spectrum extended 
morphology that is incompatible with a multiply lensed system. Thus, while these observations do 
not rule out microlensing, no macrolensing is observed. 

Subject headings: BL Lacertae Objects — AGN — Gravitational Lenses 



1. INTRODUCTION 

BL Lacertae objects are an extreme type of Active 
Galactic Nuclei (AGN), whose hallmark is their "fea- 
tureless" optical spectra (rest frame Wx < 5A for any 
emission lines present). They are a member of the 
blazar class of AGN. And like other blazars, their ob- 
served properties are thought to due to bulk relativisitic 
outflows. It is believe d that most BL Lacs are low- 
luminosity, iFanaroff fc Rilcv ( 197^ class 1 radio galax- 
ies (FR-ls) whose jet axes arc oriented roughly towards 
the observer |Blandford & Rccs 191^. Several pieces 
of evidence support this hypothesis; e.g., most BL Lacs 
have giant elliptica l host galaxies in rich galaxy regions 
ijWurtz et alJl997i) consistent with FR-1 hosts, and most 
have a highly core-dominated radio structure that is con- 
sistent with a b eamed source (tAntonucci fc Ulvestad'' 
IT985t IPerlmmfc Stockail99l lUrrv fc Padov ani 1995). 
However, there is substantial evidence that not all BL 
Lacs follow this simple picture, especially within the 
radio-selected "1 Jansky" (IJy) sample of BL Lacs 
iRector & Stocke 2001; Stickel et al. 1991). In fact, the 
extended radio power levels for fully one-half of the BL 
Lacs in the IJy sample are too luminous to be FR-ls. 
These objects also have broad, moderately luminous op- 
tical emission lines (log L = 41 — 43 erg s~^), which are 
typically seen in FR-2s but not in FR-ls. Why FR-2- 
like sources appear in abundance in t he IJy BL Lac sam- 
ple but not in X-ray selected samples iRector et alJ2000^1 
is as yet unknown, althouEfh lR.ector fc Stockel 1)20011) dis- 
cuss several possible factors. 

One of th ese factors might be gravitational lensing. 
iStocke fc R,ecto3 l|1997|) find an excess of Mgll systems 
along the sightlines of IJy BL Lacs. Ten absorption 
systems were found in the spectra of 37 BL Lacs; and 
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five of them are very strong {Wx > lA), which is four 
to five times greater than the number expected based 
upon quasar sightlines, i.e., 2.5 to 3a gre ater than the 
expectation value from iSteidel fc SargentI l|l992() . This 
excess suggests a correlation between the presence of 
absorbing gas along the BL Lac sightline and the rel- 
atively featureless spectra of BL Lac objects, as com- 
pared to quasars. Such a correlation can be caused by 
gravitational microlensing of an AGN by stars associated 
with the absorbing gas, presumably within an in terven- 
ing foreground galaxy. lOstriker fc Vietril l)1990(l postu- 
lated that BL Lacs may simply be microlensed quasars. 
In their model, a background quasar is microlensed by 
stars within a foreground giant elliptical galaxy. The mi- 
crolensing preferentially amplifies the compact nonther- 
mal continuum over the extended emission lines regions, 
reducing the observed Wx of the emission lines. If its con- 
tinuum is sufficiently boosted by microlensing, a quasar 
may appear t o be a BL Lac obiect. While it is now 
clear that the lOstriker fc Vietril l|1990fl model does not 
explain all BL Lac objects, it may explain some. And in- 
deed, evidence of lensing exists for several BL Lac s, e.g., 
AO 235-H64 (Stickel. Fried & Kiihr 1988; Webb et alJ 
22P(D, MS 0205.7+3509 ( Stocke. Wurtz & PerlmaaJ 
1995.1 . PKS 0537-441 (Romero. Suroi & Vucetich 1995^), 
PKS 1413-H3 5 (Pcrlman et al. 1996), RGB 1745-^398 
()Nilsson et al.l Il999,l and three ot her BL Lacs in an 
HST snapshot survey of BL Lacs ijScarpa et alJ 11999). 
There is also evidence against the lensing hypothe- 
s is for some of these objects , e.g., PKS 0537-441 
^Falomo. Melnick fc Tan^ Il99l and MS 0205.7-^3509 
('Falo mo et al.l 1199'^ Further, iNaravan fc Schneid'eJ 
(1990) discuss several of these sources in a microlens- 
ing context where the foreground galaxy has a gravita- 
tional potential too shallow to produce multiple images 
via macrolensing. Thus it is not known how many BL 
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Lacs are lensed, nor how many AGN are classified as BL 
Lacs as a result of lensing. 

Unfortunately, it is difficult to test for microlensing 
on a case-by-case basis. An intervening Mgll absorption 
system is suggestive, but not conclusive, of lensing. Simi- 
larly, optical imaging can, and occasionally does, discover 
an intervening foregro und galaxy which is offset from the 
BL Lac nucleus, e.g.. iHeidt et alJ l)1998|) : although this 
too is not conclusive evidence. Further, the character- 
istic light curve of a single microlensing event is likely 
to be blended into multiple events, thereby producing 
a variability signature which is not easily differentiated 
from an intrinsically variable source, e.g. iTakalo et al.l 
()1998j) : although achromatic variability is a l ensing char- 
acteristic which may be used to differentiate (jWebb et al.l 
I2OOO1) . While it is not conclusive evidence of microlens- 
ing, the detection of macrolensing (i.e., multiple images) 
remains the only means to confirm that an object is in- 
deed gravitationally lensed. So, despite the suggestive 
evidence alone, currently there is only one confirmed 
case of a lensed BL Lac: the "smallest Einstein ring" 
source B 0218-1-357. The ring is very compact; it is un- 
resolved at arcsecond-scale resolutions, but is clearly re- 
solved wi th the VLA A-arra y at 8.4 GHz and higher fre- 
quencies l)0'Dea et al.lll99^ . This compact radio struc- 
ture is typical of gravitational lenses; e.g., all the lenses 
currently known in the Jodrell Bank- VLA astronomical 
survey have angul ar separations on the order of a few 
arcseconds or less l|King et al.lll997j) . 

Here we present a sensitive, high-resolution imag- 
ing study of all of the gravitational lens candidates 
in the IJy BL Lac sample which can be reached 
by the VLA {S > -20°). Nearly all of the BL 
Lacs in the IJy sample have now been mapped at 
high sensitivity with the VLA A-arr ay at 1.4 GHz 
(jCassaro et al . 1999; Murphv. Browne fc P erlevI Il993t 
IPerlman fc' Sto cke 1994$ iRector fc Stockell200H) . Based 
upon these maps lensing candidates were chosen with 
the following criteria. First, sources were chosen that 
have an unusual radio morphology which is suggestive 
of lensing, i.e., very distorted and very luminous ex- 
tended radio structure. The B L Lacs S4 0814-f-425 
( Murphv. Browne fc PerlevI IT993|) and PKS 2131-021 
( Rector fc Stockd 1200 ID are examples of objects with 
compact (< 10") radio structure that, if interpreted 
as classical "triple" sources, appear to have edge- 
brightened, FR-2 like lobes but with a compact, wide- 
angle tail morphology that is very unusual for FR- 
2s. Alternatively, the morphology of both sources could 
be explained as Einstein rings that are marginally re- 
solved. Secondly, sources were chosen that are unre- 
solved or marginally resolved (i.e., highly core domi- 
nated). Thirdly, we chose objects whose optical spectrum 
contains a strong Mgll absorption system. This latter 
criterion selects only high-redshift objects because only 
MgH absorption systems with redshifts z > 0.4 can be 
detected at optical wavelengths. Naturally, B 0218-1-357 
was excluded from the sample since it has already been 
established to be a lensed system. Seven other objects 
from the IJy BL Lac sample meet one or more of these 
criteria and were observed in this study (see Table ^ . 
The columns are: [1] the object name; [2-3] the emis- 
sion and absorption line redshifts; [4] the 20cm extended 
radio power (in W Hz^^); and [5] the 20cm radio core 



dominance, / = Pr.nrpJ Pp.xt- Values were obtained from 
IRector fc Stockd (j200lD and references within. 

2. OBSERVATION AND REDUCTION 

We chose to observe with the A-array with a 50 MHz 
bandwidth to maximize resolution while maintaining 
good sensitivity to extended, steep-spectrum structure. 
Each object was observed at two frequencies so that 
spectral-index maps could be generated. Spectral index 
maps will show whether an object is a typical core-jet 
source, in which case the steep-spectrum jets and lobes 
will be easily differentiated from the flat-spectrum core, 
or if it is a lensed object with multiple images with 
the same spectral index. Although not attempted here, 
polarization-sensitive observations can provide an addi- 
tional test, as the polarization vectors of a lensed image 
should be unaffected by the gravitational field of the lens 
(|Dver fc ShavedlTgg^ . Thus, the lensed images of each 
source component should also have the same polarization 
percentage and position angle, once corrected for Fara- 
day rotation due to the different light paths. One caveat 
for both spectral-index and polarization measurements 
is that BL Lacs are known to vary both quantities on 
timescales as short as days. Thus, in principle multi- 
ple images of an object could display different spectral 
indices and polarizations due to the difference in light 
travel times if the variability is intrinsic and not due to 
the lens. 

Lensing candidates that are resolved at 1.4 GHz and 
have a largest angular size greater than 5" were observed 
at 4.86 GHz (C-band) and 8.46 GHz (X-band) to re- 
solve individual components while maintaining excellent 
sensitivity to faint extended structure. Smaller sources 
were observed at 8.46 GHz and 14.94 GHz (U-band) to 
improve the likelihood of resolving discrete components. 
The resolution of these maps are ten times greater than 
maps previously available for these sources. 

The seven lensing candidates from the IJy BL Lac sam- 
ple were observed with the VLA A-array on 20 Septem- 
ber 1999 and 24 September 1999. Four or five scans of 
4-8 minute durations were made for each source at each 
frequency. Scans were spaced to optimize coverage in 
the {u, v) plane. The resulting maps have sensitivities of 
roughly 0.1 to 1 mjy beam~^. 

Epoch 1995.2 VLA values were used to flux calibrate 
the maps using multiple observations of 3C 286. Since 
these sources are highly core-dominated, a point source 
model was assumed to start the self-calibration process. 
Phase-only self-calibration in decreasing solution time in- 
tervals was used for the first four iterations. Amplitude 
and phase self-calibration were then used until the maxi- 
mum dynamic range was achieved, usually requiring only 
one or two more iterations. The AIPS task IMAGR was 
used to generate the maps and clean components. Ro- 
bust weighting (ROBUST = 0.5) was used to achieve a 
smaller beam FWHM with onl y a 10-12% i ncrease in 
noise over natural weighting; see lBrigg^ l)19 9^ for an ex- 
planation. The core flux densities were measured by fit- 
ting the core with a single Gaussian with the synthesized 
beam's parameters. The extended flux was determined 
by measuring the total flux density with a box enclos- 
ing the entire source and then subtracting the core flux 
density. Throughout this paper we assume S oc 

In addition to the radio observations described above. 
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Hubble Space Telescope and ground-based optical data in 
the literature l|Pursinio et al.ll2002t IScarpa et alJl200(t 

were searched for optical counterparts to bright spots in 
the radio structure. However, no counterparts to the 
structure observed in the resolved sources were detected. 
If these are indeed lensed sources, optical counterparts 
should be present and be sufficiently bright to distort 
the observed the optical radial profiles of these sources. 
However, in all cases the objects are unresolved at optical 
wavelengths. 

2.1. Summary of Radio Properties 

The VLA continuum and spectral index maps are 
shown in Figures ^ through |H1 A summary of the radio 
properties is given in Table |21 The columns are: [1] the 
object name; [2-3] the 4.86 GHz VLA core and extended 
flux densities (mJy); [4-5] the 8.46 GHz VLA core and 
extended flux densities (mJy); and [6-7] the 14.94 GHz 
VLA core and extended flux densities (mJy). The er- 
ror estimates for extended flux densities are based upon 
the solid angular extent of the source. For unresolved 
sources, an upper limit on extended flux is conservatively 
estimated to be lOcr because faint extended flux may be 
spreaded over many beams. 

2.2. Discussion of Individual Sources 

S5 0454-1-844: This object is unresolved in 8.46 GHz 
and 14.94 GHz maps to the 0.1 and 0.4 mJy beam~^ 
level respectively. It is also unresolv ed in a 1.4 GHz, 
A-array map to the 0.1 mJy beam~^ f Uector ^, Stocka 
HLQIJ. The unresolved nature of this source is not unex- 
pected due to a redshift lower limit of z = 1.340 that is 
based upon a Mgll absorption system ijStocke fc RectoJ 
Iil97). 

PKS 0735-1-178: This object h as a in tervening ab- 
sorption system at z = 0.424 (jCarswe ll et al. 1974; 
lR,ector fc Stockgi200l . The 8.46 GHz and spectral in- 
dex maps are overlaid in Figure^ The 14.94 GHz map is 
shown in Figure [3 The 8.46 GHz map shows an unusual 
jet-like morphology to the south and west. The spectral 
index map confirms that PKS 0735-1-178 consists of a 
flat-spectrum core (af4'g4 = —0.24) and steep-spectrum 
(q^i4*94 ^ —1) jets, not multiple images of the core. 

S4 0814-H425: The 1.4 GHz map of this object in 
iMurnhv. Browne fc Perle\^ l)1993t) is highly suggestive of 
a lensed system, with possible multiple images of the core 
to the northeast and northwest. However, 4.86 GHz and 
8.46 GHz maps of this object (Figure 13 and Figure 0} 
show a jet which originates north of the core and arcs to 
the northeastern lobe. The spectral index map (Figure|2l 
confirms that the source consists of a flat-spectrum core 
(a|;4g = +0.1) and steep-spectrum (afilg < —1) jets. 
The morphology of this source is very unusual; however 
projection effects may be causing an intrinsically straight 
jet appear to bend dramatically. 

PKS 0823-1-033: This object is unresolved in both 
the 8.46 GHz and 14.94 GHz maps to the 0.3 and 
0.8 mJy beam~^ level respectively. The object is 
also unresolved in the 1.4 G Hz, A-array map in 
iMurphv. Browne fc PerlevI l)1993() . 

PKS 1749-1-096: This object is unresolved in both 
the 8.46 GHz and 14.94 GHz maps to the 0.8 and 0.9 
mJy beam~^ level respectively. It is also unresolved 
in a 1.4 GHz, A-array map to the 0.1 mJy beam~^ 



ijRector fc Stockdl200ll) . The modest tentative redshift 
(z = 0.320:) and the fact that the host galaxy is eas- 
ily observable ijWurtz et al.lll997|) make its very high ra- 
dio core dominance very unusual. However, the obser- 
vational limits allow that this source could be a typical 
FR-1 with a somewhat lower extended radio power level 
(log P20cm< 24.5 WHz-i). 

S4 1749-1-70 1: This ob ject was shghtly resolved in 1.4 
GHz map (Re ctor fc Sto ckc 200l]). The 8.46 GHz map 
(Figure|3J) shows a halo morphology around the core with 
a jet-like extension to the southwest. The halo is mostly 
resolved out of the 14.94 GHz map (Figure ^ however 
the extension to the southwest is still visible. The spec- 
tral index map (Figure ISJ shows a a flat-spectrum core 
(q^i4*94 — —0.2) and steep-spectrum (ai4*94 < —1) halo, 
indicating that this halo is not an Einstein ring image of 
the core. 

PKS 2131-021 : The 1.4 GHz map of this object in 
iRector fc Stockd l|200H ) is highly suggestive of a lensed 
system, with possible multiple images of the core to the 
east and southeast. However, 4.86 GHz and 8.46 GHz 
maps of this object (Figures [3 and ISJ show a very in- 
teresting morphology, with a jet that originates to the 
south-southeast of the core. The jet appears to be sig- 
nificantly distorted and may be helical. There also ap- 
pears to be lobes roughly 5" to the east and southeast 
of the core. It is not clear how these jets propagate to 
the lobes, although projection effects may cause a single 
lobe to appear as the two seen. But significant bending is 
still required, even with the severe projection effects ex- 
pected in BL Lacs, to cause both lobes to appear on the 
same side of the core. The spectral index map (Figure!?)) 
shows no evidence of multiple images of the core; and 
the jet and lobes have indices (a| < —1) much steeper 
than the core (a|:|^ = +0.0Q). A MERLIN-f EVN map of 
this object shows a jet component 50 mas east-southeast 
from the core, and multiple components extending to the 
south o ut to 350 mas from t he core, suggesting a helical 
motion l|Cassaro et al.ll2002|) . 

3. DISCUSSION AND CONCLUSIONS 

We have presented high-resolution radio continuum 
and spectral index maps for seven gravitational lens- 
ing candidates from the IJy BL Lac sample. All of 
the resolved sources show a flat-spectrum core and 
steep-spectrum extended radio morphology, as expected 
for core-jet sources. We note that S4 0814-1-425 and 
PKS 2131-021 have very similar, and very unusual, radio 
morphologies in that both appear to have radio lobes on 
the same side of the core. This could be caused by an 
intrinsically bent source seen close to the line of sight, al- 
though the required geometry would be rare. If so, these 
sources should have very high apparent luminosities. 
This is certainly the case for PKS 2131-021, which has 
the highest core and extended radio powers of any BL Lac 
object in the IJy samp le, discounting uncertain values 
llR,ector fc Stockdl200lli . The redshift for S4 0814+425 is 
currently unknown, although it must be z > 0.5 because 
it is unresolved at op tical wavelengths ijPursimo et alJ 
120021 lUrrv et. allll99 9*) ; thus it too must have a high ap- 
parent radio luminosity. VLBI maps of these sources 
indicate they are highly core-dominate d, supporting 
the hypothesis they are highly beamed llCassaro et alJ 
I200I iGabuzda. Pushkarev fc Cawthornd 12000(1 . Unfor- 
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tunately, further comparison is impossible without a firm 
redshift for S4 0814+425. 

None of the maps show evidence of an Einstein 
ring or multiply-imaged components; however it is im- 
portant to note that a foreground galaxy many not 
cause macrolensing if it has a low surface mass density 
ijNaravan fc Schneided riOgO). Thus, it is possible that 
a foreground galaxy may cause significant microlensing, 
brightening the continuum relative to the emission-line 
regions, without associated macrolensing. The converse 



is also true ijAbraham et al.lll993tl . Thus, these obser- 
vations do not rule out the possibility of lensing; but 
they argue strongly against macrolensing. The resolved 
sources show highly distorted jets and misaligned lobes. 
These unusual morphologies support the hypothesis that 
these BL Lacs are highly beamed sources, with radio 
morphologies that are highly distorted due to projection 
effects. 
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Object Zem Zabs log Pext f 



S5 0454+844 • • • 1.340 > 27.31 > 55.5 

PKS 0735+178 • • • 0.424 > 27.30 231 

S4 0814+425 • • • > 116.5 

PKS 0823+033 0.506: • • • 27.31 151.7 

PKS 1749+096 0.320: • • ■ 26.59 > 145 

S4 1749+701 0.770: • • ■ 27.23 15.8 

PKS 2131-021 1.285 • • • 27.99 6.9 
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Table 2. Observed Radio Properties of the Sample 



Object 


4.86 GHz 


8.46 GHz 


14.94 GHz 


/core (mjy) fext (mjy) 


fcore (mJy) 


fext (mJy) 


fcore (mJy) 


fext (mJy) 


S5 0454+844 




233.3 


< 0.5 


216.5 


< 1.5 


PKS 0735+178 




975.5 


32.2 + 3.3 


838.9 


29.3 + 6.6 


S4 0814+425 


1011.9 33.0 + 4.2 


1072.3 


18.4 + 6.7 






PKS 0823+033 




1452.9 


< 1.4 


1547.6 


< 1.8 


PKS 1749+096 




3794.6 


< 2.0 


3622.3 


< 2.7 


S4 1749+701 




520.9 


31.0 + 1.1 


459.9 


18.6 + 5.8 


PKS 2131-021 


1732.1 70.8 + 8.0 


1802.0 


36.3 + 6.7 








Fig. 1. — VLA 8.46 GHz map (contour) and spectral index map (grayscale) of PKS 0735+178. The beam is shown in the lower left 
corner. The contour levels arc 0.02, 0.05, 0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 50 and 100% the peak flux of 9.755 x lO'^ Jy beam"!. The spectral 
index range is from —1.5 to 0. 
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Fig. 2.— VLA 14.94 GHz map of PKS 0735+178. The beam is shown in the lower left corner. The contour levels are 0.05, 0.1, 0.2, 0.5, 
1, 2, 5, 10, 20, 50 and 100% the peak flux of 8.389 x lO"! Jy beam-i. 
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Fig. 3. — VLA 4.86 GHz map (contour) and spectral index map (grayscale) of S4 0814+425. The beam is shown in the lower left corner. 
The contour levels are 0.02, 0.05, 0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 50 and 100% the peak flux of 1.0119 Jy beam-^. The spectral index range is 
from —1.5 to 0. 
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Fig. 4.— VLA 8.46 GHz map of S4 0814+425. The beam is shown in the lower left corner. The contour levels are 0.02, 0.05, 0.1, 0.2, 
0.5, 1, 2, 5, 10, 20, 50 and 100% the peak flux of 1.0723 Jy beam"!. 




Fig. 5. VLA 8.46 GHz map (contour) and spectral index map (grayscale) of S4 1749+701. The beam is shown in the lower left corner. 
The contour levels arc 0.02, 0.05, 0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 50 and 100% the peak flux of 5.209 x W'^ Jy beam"!. The spectral index 
range is from —1.5 to 0. 
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Fig. 6.— VLA 14.94 GHz map of S4 1749+701. The beam is shown in the lower left corner. The contour levels are 0.02, 0.05, 0.1, 0.2, 
0.5, 1, 2, 5, 10, 20, 50 and 100% the peak flux of 4.599 x lO"! Jy beam-i. 




Fig. 7. — VLA 4.86 GHz map (contour) and spectral index map (grayscale) of PKS 2131-021. The beam is shown in the lower left 
corner. The contour levels are 0.02, 0.05, 0.1, 0.2, 0.5, 1, 2, 5, 10, 20, 50 and 100% the peak flux of 1.7321 Jy beam"!. The spectral index 
range is from —1.5 to 0. 
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Fig. 8. — VLA 8.46 GHz map of PKS 2131-021. The beam is shown in the lower left corner. The contour levels are 0.02, 0.05, 0.1, 0.2, 
0.5, 1, 2, 5, 10, 20, 50 and 100% the peak flux of 1.8020 Jy beam"!. 



